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Abstract: The AMS detector, to be installed on the International Space Station, includes a Ring Imaging 
Cerenkov detector with two different radiators, silica aerogel (n=1.05) and sodium fluoride (n=1.334). 
This detector is designed to provide very precise measurements of velocity and electric charge in a wide 
range of cosmic nuclei energies and atomic numbers. The detector geometry, in particular the presence 
of a reflector for acceptance purposes, leads to complex Cerenkov patterns detected in a pixelized pho- 
tomultiplier matrix. The results of different reconstruction methods applied to test beam data as well 
as to simulated samples are presented. To ensure nominal performances throughout the flight, several 
detector parameters have to be carefully monitored. The algorithms developed to fulfill these require- 
ments are presented. The velocity and charge measurements provided by the RICH detector endow the 
AMS spectrometer with precise particle identification capabilities in a wide energy range. The expected 
performances on light isotope separation are discussed. 



The AMS and RICH detectors 

The Alpha Magnetic Spectrometer [1] (AMS) is 
a high energy physics experiment that will be in- 
stalled on the International Space Station (ISS) by 
the year 2009, where it will operate for a period of 
at least three years. It is a large acceptance (~ 0.5 
m 2 sr) superconducting magnetic spectrometer able 
to detect in a wide kinematic range (from a few 
hundred MeV up to TeV region) singly charged 
particles, charged nuclei and 7 rays. The long 
time exposure in space will allow AMS to collect 
an unprecedented large data sample and to extend 
by orders of magnitude the sensitivity reached by 
previous experiments on dark matter and antimat- 
ter searches. In addition, the measurement of the 



cosmic -ray abundances up to the TV region and 
in a wide charge range (up to around iron) will 
contribute to a better description of cosmic ray 
production, acceleration and propagation mecha- 
nisms, essential for a full understanding of the 
background spectra on dark matter searches. Infor- 
mation about the density of the interstellar medium 
traversed by the cosmic rays and their confinement 
time can be derived from the isotopic composi- 
tion of secondary cosmic rays, produced by frag- 
mentation during the cosmic ray transport in the 
galaxy. For instance, the relative abundances of 
deuterium and helium-3 isotopes reflect the trans- 
port history along the galaxy of protons and heli- 
ums, while the beryllium- 10 radionuclide accounts 
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for the time confinement. Current measurements 
are performed at relatively low energies (T < 1 
GeV/n) and based on small statistics. 

Particle identification with AMS-02 relies on a 
very precise determination of the magnetic rigid- 
ity, energy, velocity and electric charge. In 
the AMS spectrometer, the momentum is ob- 
tained from the information provided by the sil- 
icon tracker with a relative accuracy of ~2% up 
to 10 GeV/c/n. Isotopic mass separation over 
a wide range of energies requires, in addition 
to an accurate momentum measurement, a veloc- 
ity determination with low relative uncertainty as 
Am/m = (Ap/p) © 7 2 (A/3//3). 

For this purpose, the AMS spectrometer includes 
a Ring Imaging Cerenkov detector (RICH) [2] op- 
erating between the time-of-flight and electromag- 
netic calorimeter detectors. It was designed to 
provide measurements of the velocity for singly 
charged particles with a relative uncertainty of 
0.1% and of the nuclei electric charge up to iron. 
Moreover, it will provide AMS with an additional 
contribution to the electron/proton separation. For 
the isotopic separation, the RICH detector will 
cover a kinetic energy region ranging from 0.5 
GeV/n up to around 10 GeV/n for A < 10. The 
RICH has a truncated conical shape with a top ra- 
dius of 60 cm, a bottom radius of 67 cm, and a 
total an expansion height of 47 cm. It covers 80% 
of the AMS magnet acceptance. It is a proximity 
focusing device with a dual solid radiator configu- 
ration on the top made of 25 mm thick aerogel and 
5 mm thick sodium fluoride (NaF) tiles, the later 
being crossed by <~ 11% of the events. The pho- 
ton detection is made with an array of multianode 
(16) shielded Hamamatsu tubes (R7600-00-M16), 
optically coupled to light guide acrylic pipes. 

Velocity and charge reconstruction 

Charged particles crossing the radiator material of 
refractive index n and with a velocity larger than 
1/n, emit photons. The aperture angle (9 C ) of the 
photons with respect to the radiating particle direc- 
tion depends on the particle velocity (3, 

cos6i c = -5- (1) 
p n 

The detected ring photon pattern reflects the geom- 
etry of the photomultiplier detection matrix and the 



interactions suffered by the emitted photons along 
their path to the readout matrix: radiator inter- 
actions (Rayleigh scattering, absorption), surface 
optical effects (reflection and refraction) and light 
guide efficiency. Consequently, from the point 
of view of the expected Cerenkov pattern, a typ- 
ical event will be composed of aligned photons, 
strongly uncorrected scattered photons and detec- 
tor noise. 

Two different approaches were implemented for 
the Cerenkov ring reconstruction. One was based 
on single hit reconstruction [3] and the other on 
a maximum likelihood method [4] . In the former 
method a value of (3 is reconstructed for every de- 
tected hit. The method is purely geometrical and 
relies in a set of analytical equations that relate the 
detection point with the Cerenkov angle (9 C ) and 
the particle coordinates. The possibility of the pho- 
ton being reflected is taken into account. Next, the 
most probable cluster of hits is searched and the 
final velocity is computed as a mean of the cluster- 
ized hit (3 values, weighted with measured signal 
amplitude (photon multiplicity). 

In the other reconstruction approach, the algo- 
rithm incorporates a probability density function 
for the detected hits. The residuals of the signal 
hits distribute according to a double gaussian func- 
tion whose widths are directly related to the pixel 
size and granularity, radiator thickness, chromatic - 
ity effects and aerogel forward scattering. The ex- 
istence of the second gaussian, accounting for a 
natural enlargement of the hit residuals from for- 
ward scattering, makes possible the description of 
the signal to larger hit distances and endows a bet- 
ter algorithm efficiency and a lower sensitivity to 
noisy hits. The signal probability density function 
S(r) associated to a hit i at a distance r from the 
hypothetical ring (9 C ) is as follows: 

S(r) =ai G 1 (<r 1 ) + a 2 G 2 (cr 2 ) 

where ct\ and a 2 are respectively ~ 3/4 and ~ 1/4 
and o\ and <r 2 , ~ 0.37 cm and ~ 1.35 cm. For dis- 
tances larger than ^2.1 cm, the detected hits are 
tagged as background and a uniform probability 
density function is associated, B{r) = jj, where 
D is an effective distance in the detector corre- 
sponding to 134 cm and b the background fraction 
estimated as 77%. The overall probability density 
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Figure 1 : Velocity reconstruction accuracy. 



function is therefore defined as: 

V{r) = (1-6) S(r) + B(r) 

The best value for the Cerenkov angle 9 C will re- 
sult from the maximization of a likelihood function 
C(9 C ), built from the product of the event hits prob- 
ability weighted by the photon multiplicity (for 

rip. e /hit > 1), 

N 

C(e c ) = Y[v^[n(9 c )} 

1=1 

Figure 1 shows the velocity accuracy as function 
of the electric charge obtained for simulated data 
samples. As expected, the larger is the number of 
radiated photons the better is the resolution with a 
systematic limit ~ 2 10~ 4 from the detector pix- 
elization. 

Photons are uniformly emitted along the particle 
path (L) and their number depends on the particle 
charge (z) and velocity (J3) and on the refractive 
index of the medium (n), 

-di xz L V-w^J (2) 

The fraction of photons in the ring pattern that are 
lost in every event depends on their topology (im- 
pact point and particle direction and velocity) and 
other factors. Therefore, charge reconstruction re- 
lies on the reconstructed Cerenkov angle and on 
measurements of the path length crossed by the 
particle, of the number of ring associated photo- 
electrons detected on the readout matrix and fi- 
nally on the evaluation of the photon detection ef- 
ficiency. The efficiency factors involve the radiator 
interactions, the photon ring acceptance including 
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Figure 2: Charge reconstruction accuracy. 

the mirror reflectivity and the light guide and pho- 
tomultiplier quantum efficiencies. Therefore, ap- 
plying the correction factors on an event-by-event 
basis, the detected signal for singly charged parti- 
cles (iVo) can be estimated and the charge of the 
incident particle is obtained according to: 

Z2 = ^oc^J- * . (3) 
No £tot AL sin 8 C 

Figure 2 shows the charge estimated accuracy ob- 
tained from simulated data samples for different 
nuclei. 

Detector monitoring during flight 

The accuracy of the velocity and charge recon- 
struction provided by the RICH relies on the pre- 
cise knowledge of the detector parameters through- 
out the AMS mission. For instance, the refrac- 
tive index of the individual aerogel pieces must 
be known at the 10~ 4 level not to compromise 
the velocity measurement. Similarly, the single 
PMT response must be known at the percent level 
not to degrade the charge determination capabili- 
ties. Extensive characterization of all optical ele- 
ments was performed prior to and during detector 
assembly [5]. However, the launch and, in particu- 
lar, the varying environmental conditions along the 
mission (e.g. temperatures ranging from -20°C to 
40°C) may affect the optical properties of the de- 
tector parts. Therefore, specific algorithms have 
been developed to monitor RICH optical parame- 
ters using flight data. MC simulations show that 
the radiator refractive index will be monitored to 
the required precision level using 1-day equiva- 
lent electrons or high energy protons and the PMT 
gains will be determined with the adequate preci- 
sion using a 1 -orbit equivalent (90 min.) proton 
data [6]. 
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Physics prospects: D/p separation 

The measurement of the deuteron/proton ratio pro- 
vides information on cosmic -ray production and 
propagation and is particularly challenging due to 
the low deuteron abundance (D/p ~ 1%). The ex- 
tremely accurate velocity measurement provided 
by the RICH (A/3//? ~ 10" 3 ) is crucial to reduce 
the background level. A full-scale simulation of 
the AMS detector was used to evaluate the capa- 
bilities of AMS-02 for mass separation. 

Pre-selection cuts using readings from different 
subdetectors of AMS-02 were applied to reduce 
the fraction of events badly reconstructed in mo- 
mentum. The set of events passing these cuts cor- 
responds to an acceptance of ~ 0.3 m 2 sr and ~ 0.2 
m 2 sr, respectively for protons and deuterons. 

The reconstruction of particle masses was then per- 
formed for events having a signal in the RICH de- 
tector. A series of event selection cuts were intro- 
duced, namely a minimum of 3 hits in the recon- 
structed Cerenkov ring, an upper limit on the num- 
ber of noisy hits and compatibility between two in- 
dependent velocity measurements. 

Results show that mass separation of particles with 
Z = 1 is feasible even if one species is orders 
of magnitude more abundant than the other. It is 
possible to obtain good estimates for the D/p ra- 
tio in the few GeV region (Fig. 3). D/p separation 
is possible up to Ekin ~ 8 GeV/nucleon. In the 
optimal region immediately above the aerogel ra- 
diation threshold (Ekin — 2.1 — 4 GeV/nucleon) 
rejection factors in the 10 3 — 10 4 region were at- 
tained (Fig. 4). The best relative mass resolutions 
for protons and deuterons are ^ 2% for both radi- 
ators in the regions above their respective thresh- 
olds. After all cuts, an acceptance of ~ 0.06 m 2 sr 
and ~ 0.04 m 2 sr was obtained, respectively for 
protons and deuterons for kinetic energies above 
3 GeV/n [7]. 

Conclusions 

The AMS experiment includes a Ring Imaging 
Cerenkov detector. Algorithms for charge and ve- 
locity reconstruction were developed, confirming 
the design capabilities of the instrument: up //? ~ 
0.1% for singly charged particles and a z ~ 0.2. In 
order to achieve such accuracy, an extensive char- 
acterization on laboratory of all optical elements 
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Figure 3: Reconstructed deuteron/proton ratio for 
the simulated sample. The dashed line shows the 
simulated ratio. 
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Figure 4: Rejection factor for D/p separation in 
aerogel events. 

was done and monitoring of the aerogel refractive 
index and photomultiplier gain will be performed 
during flight. Analysis based on simulated data 
samples, show that despite the large background, 
separation of deuterons from protons is possible up 
to kinetic energies around 8 GeV/n. 
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